Upper critical dimension, dynamic exponent and scaling functions in the 
mode-coupling theory for the Kardar-Parisi-Zhang equation 
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We study the mode-coupling approximation for the KPZ equation in the strong coupling regime. 
By constructing an ansatz consistent with the asymptotic forms of the correlation and response 
functions we determine the upper critical dimension d c — 4, and the expansion z — 2 — (d — 4)/4 + 
0((4 — d) 2 ) around d c . We find the exact z = 3/2 value in d = 1, and estimate the values z ~ 1.62, 
z ~ 1.78, in d — 2,3. The result d c — 4 and the expansion around d c are very robust and can be 
derived just from a mild assumption on the relative scale on which the response and correlation 
functions vary as z approaches 2. 

PACS numbers: 05.40.-a, 64.60. Ht, 05.70Ln, 68. 35. Fx 

The Kardar-Parisi-Zhang (KPZ) @ equation is a sim- 
ple nonlinear Langevin equation, proposed in 1986 as 
a coarse-grained description of a variety of growth pro- 
cesses. Due to its connection to many other important 
physical problems (such as the Burgers equation di- 
rected polymers in a random medium |^-^|, etc.), it has 
attracted much attention. When applied to a growing 
interface, described by a single valued height function 
/i(x, t) on a tf-dimensional substrate x the KPZ equation 
is: 



T){x.,t), 



(1) 



where the first term represents the surface tension forces 
which tend to smooth the interface, the second describes 
the non-linear growth locally normal to the surface, and 
the last is a noise term intended to mimic the stochas- 
tic nature of the growth process [EL We choose the 
noise to be Gaussian, with zero mean and second mo- 
ment (?7(x,t)7j(x / ,t / )) = 2DS d (x-x.')S(t-t l ). The steady 
state interface profile is usually described in terms of the 
roughness: w = (ft. 2 (x, i)) — (h(x,t)) 2 which for a system 
of size L behaves like L x f(t/L z ), where f(x) — > const as 
x — > oo and f{x) ~ x x l z as x — > 0, so that w grows with 
time like t x l z until it saturate to L x when t ~ L z . x an d 
z are the roughness and dynamic exponent respectively. 

For d > 2 there are two distinct regimes, separated by 
a critical value of the nonlinearity coefficient. In the weak 
coupling regime (A < A c ) the non-linear term is irrelevant 
and the behavior is governed by the A = fixed point. 
The linear Edward- Wilkinson equation is recovered, for 
which the exponents are known exactly x = (2 — d)/2 
and z = 2. The more challenging strong coupling regime 
(A > A c ), where the non- linear term is relevant is char- 
acterised by anomalous exponents, which are the subject 
of this paper. From the Galilean invariance Q (the in- 
variance of equation (Q) under an infinitesimal tilting of 
the surface) one can derive the relation x + z — 2, which 
leaves just one independent exponent. For the special 



case d = 1, the existence of a fluctuation-dissipation the- 
orem leads to the exact result x= 1/2, 2 = 3/2. 

While we have a satisfactory understanding of the 
KPZ equation in d — 1 Jt]|| and on the Bethe lattice [Q, 
its behavior in the strong coupling regime when d > 1 is 
controversial. Such results that are known derive mainly 
from numerical simulations || of discrete models which 
belong to the KPZ universality class ]lo| , and on several 
promising but far from conclusive analytical approaches 
JlT] [Ls[ . Most efforts are oriented towards the determi- 
nation of z as a function of d for d > 1, and to under- 
standing whether there exists an upper critical dimension 
d c above which z = 2: some analytical arguments sug- 
gest the existence of a finite d c |l4]-19[|, while numerical 
studies || and real space methods |20| find no evidence 
of a finite upper critical dimension. It is our contention 
that d c = 4, and that the numerical simulations all fail 
to see evidence for an upper critical dimension because 
in high dimensions, (d > 4), they have not been run for 
sufficiently long time for the roughness w to be in the 
scaling regime. 

We take a self-consistent approach, the mode-coupling 
(MC) approximation || , where in the diagrammatic ex- 
pansion for the correlation and response function only 
diagrams which do not renormalize the three point ver- 
tex A are retained. The MC approximation has been 
remarkably successful; it reproduces the exact values of 
the exponents in d = 1 and gives an accurate value for 
z in d = 2. Furthermore the MC equations have been 
shown to arise from the large iV-limit of an generalised 
iV-componcnt KPZ which allows in principle a sys- 
tematic approach to the theory beyond MC, by expand- 
ing in 1/N. 

Because of its relative simplicity and its seeming accu- 
racy, there have been many studies of the mode-coupling 
approach to KPZ: Hwa, Frey and Tauber [|| studied the 
d = 1 case, where it gives the exact result z = 3/2, and 
solved numerically for the scaling functions. Bouchaud 
and Cates 112] gave an approximate analytical solution 



1 



assuming simple exponential relaxation for each mode. 
Doherty et al. (l^] used instead an ansatz based on the 
form of the scaling functions in d = 1. Both [ ^2|Jl3| 
found a finite upper critical dimension, with the same 
value d c ~ 3.6. Moore et al. (l4| found an explicit so- 
lution of the MC equations for d > d c = 4 with z = 2. 
(Calculations outside the MC framework which also give 
d c = 4 include a functional RG calculation |jg| ], renor- 
malization group arguments [jl7 18 1, and an expansion 
around z — 2 |l6|]. None of these give such compelling 
arguments for d c = 4 to have won general acceptance 
either for that value for the upper critical dimension or 
even for the existence of an upper critical dimension) . Tu 
ifHf attempted the formidable task of numerically solv- 
ing the mode coupling equations by direct integration, 
and found no evidence for an upper critical dimension. 
We shall comment later on why with his particular nu- 
merical technique, problems in locating the upper critical 
dimension would arise. We believe that one of the rea- 
sons why there is a reluctance to accept that there is 



an upper critical dimension is that earlier approximate 
solutions of the MC equations yielded unconvincing non- 
integer values for the upper critical dimension (such as 
3.6...) and that the approach which did give an integer 
value of four for d c , did not, however, provide any 
insights as to the form of the solution in the physically 
accessible regime, d < 4. In this paper we shall give an 
(approximate) analytical treatment of the MC equations 
which predicts that d c is four, gives z correct to 0(4 — d) 
and provides some insights into the form of the scaling 
functions. 

The correlation and response function are defined in 
Fourier space by 



c(k,w) = <fc(k s «)/i*(k,w)>, 



G(k,w) = <$-*(]£ -k')<f" 



■w)( 



dh(k,uj) 
drjQs. ,uj ) 



where (■) indicate an average over rj. In the mode cou- 
pling approximation, the correlation and response func- 
tions are the solutions of two coupled equations, 



G-^k, u) = Gq x (k, w) + A 2 J fl / A [q ■ (k - q)] [q ■ k] G(k - q, u - Sl)C(q, Si) , 



(2ir) d 
1 2 [-d.il f d d q 



G(k, W ) = G (k, w) + \ | G(k, u>) | 2 / f J [q ■ (k - G(k - q, w - Sl)C(q, SI) , 



(2) 
(3) 



where Go(k, ui) = (vk 2 — iuj) 1 is the bare response function, and Go(k, uS) — 2D \ G(k, oS) | 2 . In the strong coupling 
limit, G(k, oj) and G(k, u>) take the following scaling forms 

G(k,w) = k-"g(u>/k*) , 
G(k,w) = k-^+ d + z ^n(uj/k z ) , 

and Eqs. (H) and (0) translate into the following coupled equations for the scaling functions n(x) and g(x): 



g 1 (x) = -ix + h(x) , 
n(x) =| g(x) | 2 I 2 (x) , 

where x — to/k z and I\ (x) and li (x) are given by 

h (af) = P Jo dO sin d - 2 9 / °° dq cos 0(cos 9 - q)q 2z ^r~ z dy g n(y) , 

h{x) = f % d9 S m d ~ 2 9f x dq(cos9~ q )*tf*-3 r -^» dyn (^) n(y), 



(4) 
(5) 



with j p = A 2 /(2 d r(^i)7r(' i + 3 )/ 2 ), r 2 = 1 + q 2 - 2qcos9. 

All the approximate analytical approaches start by 
making an ansatz on the form of the scaling functions 
n(x) and g{x). The relation z = z(d) is then obtained 
requiring consistency of Eqs. (|4|) and (||) on matching 
both sides at an arbitrarily chosen value of x. Due to the 
non-locality of Eqs. (^) and (||), the matching condition 
depends on the form of the functions n and g for all x, so 
the ansatz need to be reliable for all x. Previous ansatzes 
were unsatisfactory, since (as already pointed out by Tu 
pl| ), they did not satisfy the large x asymptotic forms 

n(x) ~ x" 1 -^, g R (x) ~ x- 1 - 2 ^, 9l (x) -> x~\ (6) 



where (3 = d + 4 — 2z and g{x) — gii(x) + igi(x). 
We will match the left and right hand sides of both Eqs. 

and Eq. (j^) as x — > oo using an ansatz that capture 
this large x asymptotic behavior. It is by doing this that 
we are able to obtain z correct to 0(4 — d). It is conve- 
nient to first write Eqs. (Q) and (||) in Fourier space: 

r%(p)=h(p), (7) 



2 (P)=-&(P), (8) 
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where I± is the Fourier transform of the real part of I\ and I2 is the Fourier transform of I2 

h (p) = 2ttP d9 sin^ 2 9 J °° dq cos 6>(cos 9 - q)q 2z - 3 gli. (pr z ) n(pq z ) , 
T 2 (p) = kPJ* d8sm d - 2 ej^° dq{cos6- q ) 2 q 2 ^ r -(d+^-2z)- n{ pq z) . 



We now make for n and g the following ansatz: 

g(p) = 2C6(p) exp(- | Dp \ 2 ' z ) , 
n(p) = Aexp(- I Bp |^/ z ) , 



(9) 
(10) 



where A,B,C,D are parameters depending on d, z and 
A. It was actually found that the direct numerical so- 
lution of the MC equations in one dimension were well- 
fitted by this ansatz ||. The same choice for g(p) (but 
not for n(p)) was also made in fl3| . The invariance of 
Eqs. (0) and (||) under the rescaling: n(p) — > n(p/a)/a 2 , 
dip) ~ ¥ d(p/ a ) (H] allows us to choose D = 1. This 
ansatz gives in x space, for large x: 

g(x) ~ const x~ l ~ 2 / z + i 2Cx~ x , 
n{x) ~ const x~ x ~^/ z , 

which have the right asymptotic properties. (This seems 
to be the key ingredient in achieving the integer value of 
4 for d c rather than a number like 3.6...). We are forced 
to set C = 1/2 to get the right large x behavior of gi. 
Then g~H is simply ie~' p ' 2/ ~. We choose to match Eqs. 
(fjj) and (||) in the limit p — > 0, which means matching 
the most divergent terms on both sides of these two equa- 
tions, and is equivalent to matching the large x behavior. 
Since | g(x) ~ 2 ~ x 2 for large x, in the small p limit, the 
left hand side in both equations (|7j |s|) i s dominated by 
the terms d 2 n/dp 2 , d 2 g~R/dp 2 , giving |22;1 : 

(2 - z)/z 2 = lim p ^o I P \ 2 - 2/z Up) , (11) 
AB^ z f3{j3 - z)/z 2 = lim p ^o | P \ 2 ~ f3/z Up) • (12) 

Performing the integrals and taking the limit one gets 
the following equations: 



PAS d d(2 - z) 1 



B 



irz 2 BI(B,z)' 



PAS d (3(/3 -z) 2 ( > 2z -W 3 



B 2 



2z-/3 



r ( a 



(13) 



(14) 



where S d = ^ dd sin d - 2 9 = 2 d - 2 T((d - l)/2)) 2 /T(d - 1) 
and I{B,z) = J °°ds(l - 2 S 2 )s 2 ^ 3 e^{~B^ z s^ - s 2 ). 
Setting d = 1, z = 3/2 in the expression for B = B(z, d) 
obtained from Eqs. ([l3]) and (jlj) we find that B ap- 
proaches 1, and then (since (3 — > 2) n(p) oc g~R.(p), 
as required by the fluctuation-dissipation theorem. In 
the limit d — > we get the correct value z — 4/3. 
As z approaches 2, BI(B,z) — > and the right 



hand side of Eq. (|T^) goes to zero. It is reasonable 
to assume that A is a finite number [^3), which forces 
B ~ (2 — z)~ 1 . From comparison with equation ( |l4| ) one 
gets 2z - f3 = 4 — d — 4(2 — z) ~ (2 - z) 2 which gives 
4 - d = 4(2 - z) + 0((2 - z) 2 ). Note that this result, not 
only fixes d c to be 4, but also provide an expansion for z 
around d~: z = 2 — (4 — d)/4. In the simplest scenario 
B = (const(2 — z)) , where the constant is fixed to 2 
by the value of B in d = 1. Eqs. ( |l3| ) and (|l4|) are then 
easily solved numerically for 1 < d < 4, giving z ~ 1.62 
in d = 2 and z ~ 1.78 in d = 3. The result in d = 2 is in 
good agreement with the values found in the literature 
from simulations (which of course do not use the MC ap- 
proximation) and also with the direct numerical solution 
of the MC equations by Tu, which we suspect may be 
accurate for d < 2. 

The results d c = 4 and^ = 2-(4-d)/4+0((4-d) 2 ) are 
very robust: assuming a simple scaling form for n(p) as 
z approaches 2, they can be derived without any specific 
assumption on the form of the scaling functions, other 
than the following weak requirements for the small p be- 
haviour of the scaling functions: g~B,{p) — C(l — kp 2 / z ), 
n(p) ~ A(l — fc'(p_B)^/ z ), where k and k' are constants 
of order 1 as z — > 2. The scaling form for n(p) is: 



lip) - Af(Bp) , 



(15) 



where B = (2(2 - z))-\ /(0) = 1, and A/(0) is finite as 
z — > 2 J24j] (/(#) denotes the inverse Fourier transform of 
f(p))- We fix the scale on which g varies to be of order 1 
p5j. When z is close to 2, g~R{p) is a slowly varying func- 
tion compared with n(p), which allows us to perform the 
integral h(p) to the leading order in (2 — z). Matching 
Eq. (R) at p — gives: 



PAS d _ d(2 - z) 2kz 
B ~ irz 2 f(0Y 



(16) 



Thus B ~ (2 — z) , consistent with our assumption. A 
similar analysis of Eq. (pf) in the limit p — > gives 



i? 2 



(17) 



which again leads to B 2 ~ (4z — d — 4) l , and thus to 
d c = 4, z = 2-(4-d)/4 + 0((4-d) 2 ). Eqs. (pj) and @ 
are of course consistent with Eqs. (|l3|) and (O) obtained 
from our original ansatz. 

Let us now briefly discuss what could have gone wrong 
in Tu's numerical analysis. Assume our ansatz (^) for 
g(p) is reliable for all p as d approaches 4. Then in d = 4 



3 



g(x) — (g(0) _1 — ix)^ 1 . Thus the imaginary part of g 
has a peak at x = t/(0) _1 . We broke the scale invariance 
by choosing D = 1, which gives <?(0) = 1 and lead to 
n(0) ~ (2 — z). Tu chose instead to break the invariance 
by fixing the value of n(0) to be some constant. Rescal- 
ing our solution to make it compatible with Tu's choice, 
we would get g(0) ~ (2 — z), which causes the peak in the 
imaginary part to run to infinity, making it very difficult 
to get an accurate evaluation of the integral I\{x). It 
would be very interesting to repeat Tu's analysis in the 
light of our findings on the scaling form of n(x) near the 
upper critical dimension. 

The validity of our results beyond the mode-coupling 
approximation is an important consideration. We expect 
that the result d c — 4 will hold beyond mode-coupling, 
as argued in Ref. ]l4| ]. However, we know of no argument 
that would imply that z ~ 2 — (4 — d)/A would remain 
true outside the MC approximation i.e. there might be 
corrections at 0(A — d) to the value of the dynamical 
exponent z from the diagrams not included in the MC 
approximation. 

To summarise, we have developed a self-consistent ap- 
proach to the strong coupling regime of the KPZ equa- 
tion. We constructed an ansatz for the scaling func- 
tions which is consistent with known asymptotic results, 
and gives naturally an integer upper critical dimension 
d c — 4. ft also provides for the first time a satisfactory 
expansion of z around d c . Solving for z = z(d) gives 
good agreement with previous studies for all 1 < d < 2. 
We showed also how the same results can be obtained 
with much weaker requirements on the form of the scal- 
ing functions, if one assumes a scaling form for n. We 
found a possible mechanism for the failure of Tu's nu- 
merical analysis as d approached d c . Our results together 
with previous numerical MC studies for d = 1 and d = 2, 
f|[2l| , and with the analysis above d = 4 of Ref. [jlJJ , give 
a semi-quantitative picture of the behavior of the strong- 
coupling regime of the KPZ equation for all dimensions 
d. 
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